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Abstract 

The Hamiltonian structures of several hybrid kinetic-fluid models are identified explic- 
itly, upon considering coUisionless Vlasov dynamics for the hot particles interacting with 
a bulk fluid. After presenting different pressure-coupling schemes for an ordinary fluid 
f~^ ' interacting with a hot gas, the paper extends the treatment to account for a fluid plasma 

^SJ , interacting with an energetic ion species. Both current-coupling and pressure-coupling 

MHD schemes are treated extensively. In particular, pressure-coupling schemes are shown 
^ , to require a transport-like term in the Vlasov kinetic equation, in order for the Hamilto- 

nian structure to be preserved. The last part of the paper is devoted to studying the more 
general case of an energetic ion species interacting with a neutralizing electron background 
(hybrid Hall-MHD). Circulation laws and Casimir functionals are presented explicitly in 
r^ ', each case. 
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1 Introduction 

1.1 Hybrid plasma models 

While fluid models are widely successful in plasma physics, kinetic eff'ects have been shown 
to be relevant in many situations involving coexistence of cold plasmas and energetic hot 
particles. The latter play an important role in different contexts, ranging from fusion research 
[36] to astrophysical plasmas [45]. The need for multiscale models that accommodate the 
statistical kinetic effects of hot particles is a subject of current research, particularly involving 
computer simulations [35]. These models are usually realized by following a hybrid philosophy 
that couples ordinary fluid models to appropriate kinetic equations governing the phase-space 
distribution of the energetic particle species. 

The two main research directions involve either the coexistence of a fluid MHD plasma 
with an energetic ion component [36] or the coexistence of an electron neutralizing background 
with an energetic ion component [45, 47] (hybrid Hall- MHD). The first direction splits into two 
possible approaches: the current-coupling scheme [36, 1, 43] and the pressure-coupling scheme 
[36, 11, 7, 42], depending on how the fluid equation is coupled to the kinetic equation for the 
hot particles. 

In particular, the pressure-coupling scheme possesses some variants that will be discussed in 
this paper. In all the variants appearing in the literature, the kinetic-fluid coupling occurs solely 
as an extra pressure term in the fluid momentum equation, while the hot particles are affected 
by the cold background plasma only through the electromagnetic field terms appearing in the 
kinetic equation. The fiuid plasma density is conventionally transported by the background 
fluid and this transport generates the usual (barotropic) pressure effects in the evolution of the 
fluid momentum. On the other hand, the kinetic pressure corresponding to the hot particles 
has never been considered as arising from fluid transport terms appearing in the accompanying 
kinetic equation. The absence of these fluid transport terms destroys certain relevant properties 
of ordinary fluid models, such as the celebrated Kelvin circulation theorem that is of paramount 
importance, for example, in geophysical fluid dynamics. Also, while the momentum equation 
usually conserves the prescribed energy, the absence of transport terms in the kinetic equation 
breaks the corresponding Hamiltonian structure, which is of interest in various plasma models 
[27, 31, 28, 3]. However, even when considering hypothetical transport terms, their inclusion in 
the kinetic equation cannot be achieved by ad hoc arguments and their insertion must follow 
from general principles. An interesting work in this direction was carried out by Fla [0], who 
considered a linearized approach. Fla's results follow from the geometric properties underlying 
kinetic Vlasov-type equations and in particular their Lie-group symmetries [25, 44]. 

The pressure-coupling scheme has been adopted in many simulation codes, probably be- 
cause of its simplicity (compared to the current-coupling scheme). On the other hand, this 
scheme assumes that the averaged kinetic momentum of the hot particles is negligible at all 
times. This assumption is not required by the current-coupling scheme, which then possesses 
a wider range of validity. The current-coupling scheme couples the fluid and kinetic equations 
through the Lorentz force that is exerted by the hot particles on the background fluid plasma. 
Although this scheme has been used in [1, 43], it remains less popular in the plasma physics 
community. Nevertheless, this paper stands in favor of this model, because of the strong 
assumption on the hot particle momentum that is required by the pressure-coupling scheme. 

All existing hybrid codes usually replace the phase-space kinetic equation for the energetic 
component by particle dynamics. The distribution function is then reconstructed in order to 
calculate the macroscopic variables, such as pressure and momentum. In particular, the hot 
particles are often advanced by appropriate gyrokinetic equations [23, 13, 2]. Nevertheless, this 
paper will describe the energetic ions in terms of the more general Vlasov kinetic equation, 
whose Hamiltonian formulation is easier then its gyrokinetic counterpart. The particle picture 



can always be recovered by the (^-like particle solution, which is well known to preserve the 
Hamiltonian structure [19]. 

1.2 Hamiltonian methods in plasma physics 

The use of Hamiltonian methods in plasma physics goes back to the early 80's, when many well 
known plasma models were found to possess a Hamiltonian structure. The latter is composed 
of a Poisson bracket and a Hamiltonian function(al), which unambiguously identifies the total 
energy of the system. For example, celebrated Hamiltonian structures are those for multi-fluid 
plasmas [41, 40], for magnetohydrodynamics [34, 17] and for the Maxwell- Vlasov equations 
[33, 25, 44]. Another important discovery was the Hamiltonian formulation of guiding-center 
motion, which led to the modern theory of gyrokinetic equations [13, 23, 2]. The development 
of Hamiltonian techniques is still an active area of research, as shown by the recent works on 
magnetic reconnection and gyrofluid models (see [31] and references therein). 

The Hamiltonian formulation is particularly advantageous when it is accompanied by a 
variational principle on the Lagrangian side, which is the case for most plasma models [31]. 
When this happens, this is due to the fact that the Poisson bracket arises naturally from 
the relabeling symmetry of the plasma [27]. The advantage of such a Hamiltonian structure 
resides in several aspects, whose most celebrated example is the possibility of carrying an 
exact nonlinear stability analysis [IS]. This analysis is made possible by the existence of 
invariant functions (Casimirs) that arise from the Poisson bracket. On the Lagrangian side, 
Noether's theorem produces conserved quantities, thereby establishing the existence of explicit 
circulation laws that are a fundamental feature of fluid models. Moreover, the variational 
principle associated to the Hamiltonian structure can be approached to apply asymptotic 
methods or Lagrangian averaging (or Lie-transform) techniques. Some recent reviews on these 
topics, with special emphasis on plasma physics, are available in [3, 31, 28]. 

In certain cases, the Hamiltonian structure does not possess an immediate variational 
principle formulation in terms of purely Lagrangian variables. For example, although the 
Maxwell- Vlasov system [33, 25, 44] possesses several variational formulations [37, 38, 39, 46, 24], 
these are based either on Eulerian variables or on a mixture of Eulerian and Lagrangian 
variables. On the other hand, purely Lagrangian variables were used in [5], upon suitably 
modifying Low's Lagrangian [24]. The method in [^)] is general enough to be applicable to any 
Vlasov-type system. This paper relies on this approach to produce the Lagrangian formulation 
of the Hamiltonian plasma models that will be presented. 

1.3 Goal of the paper 

This paper aims to present the explicit Hamiltonian structure of several hybrid kinetic-fluid 
models, for either ordinary barotropic fluids, MHD and Hall MHD. In particular, this paper de- 
rives a whole class of hybrid kinetic-fluid models by making use of well established Hamiltonian 
Poisson bracket methods. Upon taking the direct sum of ordinary Poisson bracket structures, 
the hybrid models are derived by simply transforming the bracket appropriately, so that the 
Hamiltonian structure is always carried along in a natural fashion. The current-coupling MHD 
scheme is found to possess an intrinsic Poisson bracket, which is presented explicitly. On the 
other hand, pressure-coupling schemes is derived by neglecting the momentum contribution of 
the hot particles directly in the form of the Hamiltonian, rather than approximating the equa- 
tions of motion in their final form. This key step of approximating the Hamiltonian instead of 
the equations leads to the preservation of all the usual properties of fluid models, such as cir- 
culation theorems and even helicity preservation. All the considerations in this paper restrict 
to consider non-collisional Vlasov dynamics and barotropic fluid flows. The more general case 
of adiabatic fluid flows accounting for specific entropy transport is a straightforward extension. 



After presenting the most basic (pressure-coupling) hybrid models for ordinary neutral fluids, 
the paper formulates hybrid models for electromagnetic fluids. Hybrid versions of magnetohy- 
drodynamics are presented, together with its associated circulation theorems. The last part of 
the paper is devoted to presenting the Hamiltonian structure of the hybrid Hall-MHD scheme 

[45, 47]. 

Remark 1 (The role of Lie symmetries) The whole paper makes use of well established 
Lie-symmetry techniques to produce the various Hamiltonian structures. Indeed, the use of 
Lie-symmetry concepts (e.g. momentum maps and group actions) is of central importance to 
the treatment, although the discussion ignores all technical points. Consequently, all the result- 
ing Poisson brackets possess an intrinsic geometric nature, which goes back to the invariance 
properties underlying MHD [li'i] and the Maxwell-Vlasov system [25, 27]. Thus, the emergence 
of Lie-Poisson systems is an immediate consequence of the treatment. Restricting to consider 
these symmetric Lie-Poisson structures is a deliberate choice of the author and the possibility 
of other forms of Poisson brackets is not taken under consideration. 

2 Elementary hybrid models 

This section presents the most basic ideas underlying hybrid models. We consider two simple 
situations: the hybrid description of a single fluid and the hybrid model for an ideal fluid 
interacting with an ensemble of hot particles. 

2.1 Hybrid formulation of ordinary fluid dynamics 

The most fundamental hybrid model already arises in the well known derivation of fluid dy- 
namics from kinetic theory [22]. Upon restricting to the non-collisional case, one considers a 
Vlasov equation of the form 

where p is the momentum coordinate (the single-particle mass is taken to be unitary for 
simplicity) and /(x, p,t) is the Vlasov distribution on phase-space. The collective force field 
F = — V(/>(/) arises from a potential 4>, which is usually a linear functional of /. Equation (1) 
conserves the total energy 

^(/) = ^//(|p|' + ./.(/)) d^xd^p (2) 

and possesses the Poisson structure [30, 31, 27, 28, 32, 29] 

where {•, •} denotes canonical Poisson bracket on phase space. Thus, equation (1) becomes 

The equations of fluid dynamics follow by closing the equations of the moments [12] 

^(x,i)= /p«^/(x,p,t)d3p (5) 



where p®^ is the s-th tensor power on the momentum coordinate. The moment hierarchy 
is conventionally closed to consider the first two moments {Aq,Ai) as the basic dynamical 



variables, while the second-order moment dynamics is usually assumed to depend on (Aq, Ai). 
Evidently, n = Aq is the particle density, while K = yli is the averaged kinetic momentum, 
that is the fluid momentum. We notice that the Hamiltonian (2) can be rewritten in terms of 
the first two moments as 



H{n,KJ) = ^ If 



K 

P 

n 
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d^x + $(n), 



(6) 



2 J n 
where $ denotes a potential energy functional depending exclusively on the particle density n. 



Remark 2 (Mean and fluctuation terms in the Vlasov Hamiltonian) Notice that the 
above form (6) of the Vlasov Hamiltonian (2) has the advantage of splitting explicitly the 
mean and fluctuation parts of the kinetic energy, corresponding to the second and first term 
respectively. When the fluctuation part is absent, the system reduces to a 'cold plasma ' system, 
which usually arises from the simple moment closure f = n5{p — K/n). On the other hand, 
in certain situations involving small amounts of energetic particles, the energy contribution of 
the averaged momentum K is often neglected. Then, the second term in (6) is small and all 
the kinetic energy is concentrated in the fluctuation part, which is the first term. As we shall 
see, this point is of central interest in hybrid models for plasma physics. 

In order to express the Vlasov equation (1) in terms of its macroscopic fluid quantities, one 
inserts (7) the chain rule formula 



6F _ 5F 5F 6F 

Jf ~ J^^'^'^^Jf 



(7) 



in the Poisson bracket (3), thereby obtaining the Poisson structure 
{F,G}(n,K,/)= /K- 



+ / 
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5F 6G 



d^x 



SF 5G 6G ^SF\ ., 
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(8) 



Here, [X,l^] = —{X ■ \/)Y + (Y ■ V)X is minus the commutator on vector fields. Upon 
dropping the bar notation for convenience, Hamilton's equations are written as: 



dn , f 6H 

dK ,. f5H\ [^5H\ ^^ 
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(10) 

(11) 



These must be accompanied by the constraints 



K- /p/dp = 0. 



n 



fdp = 0. 



which enforce the moments (K, n) to depend on the distribution function / at all times. Indeed, 
although the hydro dynamical quantities (K, n) are mutually independent, they both depend 
on the statistical distribution /. (In more geometric terms, the above constraints identify a 
zero-level set of a momentum map, known as plasma-to- fluid map [27]). 



At this point, the Hamiltonian (6) yields 
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thereby producing the pressure term as follows: 

Proposition 1 With the notation above, the following relation holds 



^/p{/,|p-Fp}d^P = div|(p-yfVd='p. 



where we have defined the velocity Y = K/n. 

Proof. The proof is given by direct verification as follows: 
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i /p {/, IP - vp} d'p = /p {/, (i IpP - p . V + i I vp) } d'p 

= div /pp / d^p - {V ■ V)K - VF • K - (div F )K + -nV\V\^ 

= div /pp / d^p - {V ■ V)K - (div F)K 

= div (n VV + /(p -V){-p-V)f d^p j - {V ■ V)K - (div V)K 

= divy"(p-F)(p-F)/d3p. ■ 
Therefore, the hybrid equations of motion read as 

'" " (12) 

ip-Vf^fd'p (13) 

(14) 

where we have defined the potential energy density (f) = 6^/6n and V = K/n. 

At this point, one usually proceeds by invoking thermodynamic principles leading to an 
equation of state. This simplifies the pressure tensor on the right hand side of the velocity 
equation and allows to discard the kinetic effects contained in the Vlasov equation for /. This 
process leads to ordinary fluid equations for inviscid barotropic fluids. However, when kinetic 
effects cannot be neglected, the Vlasov kinetic equation must be retained in the above system, 
which then constitutes the hybrid formulation of ordinary fluid dynamics. 

We notice that the terms on the left-hand side of the kinetic equation in (11) play a key role 
in the Hamiltonian structure of the system by coupling fluid and kinetic terms in the Poisson 
bracket. Even more importantly, the term p • 5H/5Ji automatically generates the pressure 
term in the velocity equation: this presence of the pressure tensor in the velocity equation is 
not allowed when the term p • 6H/6Ji is absent in equation (11). Moreover, the same term 
ensures the validity of the Kelvin circulation theorem 

-^ / F-dx=-/ (- div [{p-V)®^ fd^p] -d^, (15) 

dt Jytiv) ht(V) \P J J 

where 7i(V) is an arbitrary loop moving with velocity V and the right-hand side evidently 
arises from the kinetic energy associated to the fluctuation velocity p — V. When this fluc- 
tuation velocity is absent in the Hamiltonian (6) the right-hand side above vanishes, thereby 
returning the usual Kelvin's theorem (f V = const. 



2.2 Hamiltonian hybrid kinetic-fluid systems 



This section considers the interaction of an ideal compressible fluid interacting with an ensemble 
of energetic particles. This basic example is of fundamental importance because it embodies all 
the main properties that will be used later to derive hybrid pressure-coupling plasma models. 
In order to derive the hybrid system, we begin by considering the following equations for 
the fluid density p and velocity u and for the Vlasov distribution /: 
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(16) 
(17) 
(18) 



Here p denotes the scalar fluid pressure while the total potential energy $(p, /) is a functional of 
both the fluid mass density and the Vlasov distribution. For instance, the case of electrostatic 
interactions involves a potential energy functional of the form [30, 27] 

^ = -\ll[lc P(x) + g^y /(x, p) d^p^ ^ l^q, p(x') + qhjfi^', p') d^p'^ d^x d^x' , 

where the charge labels c and h denote the fluid (cold) component and the energetic (hot) 
component respectively (again, we consider unitary particle masses for simplicity). Then, the 
potential energy represents the only coupling term of the above equations, which otherwise are 
completely decoupled. 

Since the Hamiltonian structures of ideal compressible fluids and kinetic Vlasov equations 
are well known [27], it is easy to see that the above system possesses a Hamiltonian formulation 
with the direct sum Poisson bracket 
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(19) 



and the Hamiltonian 



H{m, P,f) = l y"^ d^x +^Jf |p|2 d^x d^p + $(p, /) + Jp U{p) d^x , (20) 

which are both written in terms of the fluid momentum xn = pu. Notice that we have 
introduced the fluid internal energy lA such that p = p^ dU /dp. 

At this point, in the search for a hybrid kinetic-fluid model, one defines the total momentum 

M = m+ / p/d^p 



5F _ 5F 6F 
If ~^"6M^Jf' 



and expresses the equations of motion in terms of the new variable. This can be done by 
replacing the chain-rule formulas 

dF _ 6F 5F _ 6F 

'5m ~ m.' "5p~"^' (5/ ^(5M 

in the Poisson bracket (19), thereby obtaining the structure 
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(21) 



Then, upon denoting K = Jp / d'^p, the total energy expression 



H(M,p,f) = -J^ Ld'^ + -Jf\p\^d^^d^p + ^pJ) + jpUip)d'^, (22) 

yields the equations of motion 

^ + div{pu) = (23) 

-^— + uVu= — div / pp/d^p- -Vp- V-^ (24) 

dt p J p 6p 

where U = M//) and u = m/p = (M — K)//9. These last equations of motion are totally 
equivalent to the system (16)- (18), that is they yet involve no approximations. 

2.2.1 First pressure-coupling scheme 

As we remarked in the preceding section, one normally considers small amounts of energetic 
particles, so that the total energy contribution of the averaged momentum is often neglected. 
In the plasma physics literature, this assumption is made by simply replacing the convective 
term u ■ Vii by U ■ VU in the velocity equation. On the other hand, this simple step in- 
exorably destroys the Hamiltonian structure of the equations of motion, which then lack an 
energy balance equation and therefore need substantial modifications. A simple possibility of 
preserving the Hamiltonian structure is provided by the assumption M = m + nV^~min 
the Hamiltonian (22). Then, the total energy 

H{M,p,f) = -J^^d'^+-Jf\p\^d'^d^p + Hp,f) + JpU{p)d^^, (26) 

yields the equations of motion 

^ + dW{pU) = (27) 

^ + [/-VC/ = --div/"pp/d3p-lvp-V^ (28) 

dt p J p 6p 

df r r ^r, df d ^^ df 

i + l/.P-t^l + P-^-^^.^-O. (29) 

While the velocity equation is exactly the same as that of some hybrid schemes frequently used 
in plasma physics [11, 36, 7] (up to Lorentz force terms that will be included in the following 
sections), the Vlasov equation carries the extra term 

... ^^. ^^ df df dU 
{f,p-U} = U-- — ^- • li- • P 
ax dp ax 

which evidently arises from the velocity shift U = u + K/p. 

Remark 3 (Transport of the Vlasov distribution) We emphasize that the term {/, p • U} 
is required by the Hamiltonian structure of the original equations and cannot arise from direct 
approximations on the equations of motion (16) -(18). The same term was found by Fid [9], 
who used more sophisticated techniques. We notice that the term U ■ V/ leads to the inter- 
pretation of a Vlasov equation that is 'transported' along the macroscopic velocity U. This 
transport property generates the circulation force term dpf ■ VC/ • p that is related to the change 
of reference u ^ U . The kinetic-fluid interaction term {/, p • U} is the main novelty of this 
model. 



The Vlasov equation may be rewritten as (?(/+{/, 1/2 ( |p + C/| —\U\ )+(5y'<I>} = 0. Notice 
that the equation dtn + div{nU) = — divK for the energetic particle density n = Jfd^p 
involves the continuous source term divK. Thus, a slowly varying density n enforces the 
averaged momentum K to possess negligible divergence. 

2.2.2 Second pressure-coupling scheme 

So far, we only made the simple assumption M ~ m in the expression of the total energy. This 
was justified by the fact that the energetic particles exist only in small amounts, so that the 
averaged momentum K is assumed not to contribute to the total energy. However, as we have 
remarked in section 2.1, under this assumption the kinetic energy of the energetic component 
may be modified to neglect its average value |Kp/(2n). Then, the total energy becomes 
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(31) 
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The same form of the above velocity equation also appears in certain hybrid schemes [42], 
although the above Vlasov equation differs substantially from those models because the present 
system accounts for energy balance, contrarily to the pressure-coupling schemes appearing in 
the physics literature. We notice that the above Vlasov equation can be written in the form 
dtf + {/, l/2( |p — F|^ + |p + C/p — IpP — |t/p ) + <5/*^} = 0, which emphasizes the various 
kinetic energy contributions. A natural consequence of this kinetic equation is that its zero-th 
moment n satisfies the simple advection relation dtn + div(n U) = 0, so that the total density 
D = p + n satisfies the equation dtD + div(D U) = 0, as it appears in [21]. 



2.2.3 Third pressure-coupling scheme 

The advection equation for the total density D arising in the last case, suggests looking at the 
equations of motion in terms of this density variable. To this purpose, one again replaces the 
chain-rule formulas 

6F _ 6F 6F _ 5F 6F _ 5F 6F 6F 



in the Poisson bracket (19), thereby obtaining the following Poisson bracket 
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which is identical to (8). Under the assumption of a rarefied energetic gas, one has D c^ n 
and each of the Hamiltonians used above yields a slightly different system. For example, upon 
replacing p by D, the energy in (30) produces the equations 

dD 

— + div(i? [/) = (34) 



^ + U-VU = -^ divy"(p -V){p-V)f d^p - ivp - V^ (35) 

|,{;,p.(„_V).i(,.P-|.P)}.pl-|;|.|^0. (se) 

which involve no substantial modification with respect to the previous hybrid schemes. 

2.3 Conservation laws 

The Hamiltonian structure of fluid models has the main advantage of preserving most of the 
geometric features possessed by ordinary ideal fluids. One of the most significant features is 
the Kelvin circulation theorem 

^ / w • dx = (37) 

that holds for the velocity equation (17). Here the curve 74(1*) is a loop moving with velocity 
u. The hybrid schemes available in the physics literature often lack this fundamental feature, 
which is instead available for the hybrid models presented above. Indeed, the momentum 
shift M = m + /p/d'^p together with its resulting Poisson bracket (21) guarantee the Kelvin 
circulation theorem 

It is important to observe that the circulation law (38) is different from (37), since the loop 7t 
moves with velocity u in (37) and with velocity U = u + p~^ fp f d^p in (38). This is due to 
the fact that modifying the fluid kinetic energy in the Hamiltonian (22) yields a different fluid 
velocity. 

Remark 4 (Relation to turbulence models) Notice that exactly the same phenomenon 
occurs in certain regularization models for fluid turbulence [ii], such as the Navier-Stokes-a 
model [10]. These models have also been called "Kelvin- filtered" turbulence models. Another 
property of the preceding pressure- coupling schemes that is shared with Kelvin- filtered mod- 
els is that the fluid vorticity u = V x u is transported along the modified velocity U, that 
is dfUJ -\-'V X (U X Lj) = 0. This equation evidently possesses vortex filament solutions that 
project to point vortices in two dimensions. The appearance of this close analogy between 
pressure- coupling hybrid schemes and fluid turbulence models looks promising for further in- 
vestigation, especially in relation to the energy splitting appearing in the Hamiltonian (6). 

The circulation conservation (38) is valid in any of the pressure schemes presented previ- 
ously. In the third pressure-coupling scheme, the Kelvin circulation theorem becomes 

(u - (d- ff d^p] /p/d^p) .dx = 0. 



dt j^^(u) 
Moreover, the previous hybrid schemes possess the helicity conservation relation 

^J Uj--Jpfd?p\ .Vx Uj--Jpfd^p\d?i, = Q. 

10 



Notice the above conserved quantity is a Casimir of the Poisson bracket (21) and in two 
dimensions, this Casimir can be used to approach the problem of nonhnear stabiHty, fohowing 
the reference [18]. Another Casimir of the Poisson bracket (21) is given by C = / ^{f) d^xd^p, 
for an arbitrary function $. This is the Casimir typically associated to Vlasov dynamics and it 
includes the total mass //d^xd^p. The search for other Casimir functionals is left for future 
investigation. 

3 Hamiltonian hybrid MHD models 

The purpose of this section is to formulate a kinetic-multifluid model that leads to the formu- 
lation of hybrid MHD schemes. As we shall see, the pressure-coupling schemes can be obtained 
by repeating systematically all the steps outlined in the previous sections, upon including elec- 
tromagnetic fields appropriately. On the other hand there is another relevant hybrid MHD 
scheme, which involves a current coupling. This section derives both current and pressure 
coupling schemes and provides their Hamiltonian structures. At the end of this section, the 
two schemes are compared in terms of their energy-balance properties. 

3.1 Kinetic-multifluid system 

The equations of motion for a multifluid plasma in the presence of an energetic component 
read as 

Ps-^ + Ps {us ■V)us = asPs (E + Us xB) - Vp^ (39) 

^ + div ip,u,) = (40) 

f + J^.|^ + ,jE + ^xB).f = (41) 

ot rrih ax \ rrih J op 

Moeo-^ = V X B - /io ^ asPsUs - Poah / P / d^p (42) 

- = -VxE (43) 

eodw'E = ^asPs + qh fd^P, divB = (44) 

s ■^ 

where a^ = qs/mg is the charge-to-mass ratio of the fluid species s, while ps and Ug are its 
mass density and velocity. Notice that in the above system we have restored the hot particle 
mass ruh- When the energetic component is absent, the Hamiltonian structure of this system 
was discovered by Spencer [40, 41] and analyzed further in [27]. Combining this Hamiltonian 
structure with that of the Maxwell- Vlasov system [25, 44, 33] yields the following Poisson 
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bracket for the kinetic-multifluid system 



f ^ I \ Am. A/1. Am. An. / 



5ms 6ps Suis 6psJ 

V- /■ [ 5F 5G 5G 5F ^ 6F 5G \ . 
+ > / asPs (^-T^-^-T^ + ^-l: — Xt — dx 
^-^ J \ova.s OLJ orris OIJ ouig ouis J 

(■ f d 5F 6G d 6G 5F ^ f d 5F d 5G\ , ,o o 



\dp6f 6B dp6f 6B \dp 6f dp 6f 

(■ (5F JG 5G JF\ ,3 ,,^, 

where evidently m^ = psUs- To complete the Hamiltonian structure of the kinetic-multifluid 
model, one writes the Hamiltonian 

i7(rrr„/,„/,D,B) = i^ A ^:^ d^x + -^ / / |p|2 d^xd^p + ^ fpsU{ps)d^^ 

^ g J Ps ■^m/i J s 

+ Ul [\nfd'^+- [\B\'d'^) , (46) 

where the electric induction D = eoE is (minus) the conjugate variable to the vector potential 
A, as explained in [14, 15]. 

Remark 5 (Extendirrg to several errergetic componerrts) The above kinetic-multifluid 
system can be easily extended to the case of several energetic components, each described by its 
own Vlasov equation. Then, each of these Vlasov equations can be rewritten in terms of the 
first two moments (n, K) in order to isolate kinetic pressure terms. A similar approach was 
followed in [8]. 

3.2 Current-coupling hybrid MHD scheme 

In usual situations, one is interested in one-fluid models. Thus, it is customary to special- 
ize to the two-fluid system and to neglect the inertia of one species (electrons). This last 
approximation is equivalent to taking the limit m2 — >■ for the second species in the total 
fluid momentum equation. Under this assumption, the sum of the equations (39) for s = 1,2 
produces the equation 

pi— — -I- pi {ui ■V)ui = {aipi + 02/92) E + {aipiui + 02/921*2) x B - Vpi (47) 

Also, upon assuming neutrality by letting eo —5- 0, the electromagnetic fields satisfy the equa- 
tions 

V asPsUs = —V xB- Oh pfd^p (48) 

^ f^o J 

f)B 

— = -V X E (49) 

^asPs = -qh I fd^p, divB = (50) 
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Then, equation (47) becomes 

p^ + p{u ■\/)u = - Lh j f d^p]B + (—V xB - ahjpfd^p] xB-pVp, (51) 

where we have dropped labels for convenience. Finally, inserting Ohm's ideal law E+ix x B = 0, 
the kinetic two-fluid system becomes 

p^ + p{u-V)u= LhU ffd^p- aft /p/d3p + — VxBJ xB-pVp (52) 

fin 

(53) 

u) X B • -^ = (54) 



+ div (pu) = 




Of ^ P df 1 




^ = V X (ti X B) 





dp 



(55) 



which is identical to the current-coupling hybrid scheme presented in [11, 36, 1], except the 
fact that particle dynamics is governed by the Vlasov equation rather than its gyrokinetic 
counterpart. Notice that the above system does not make any assumption on the form of the 
Vlasov distribution for the energetic particles. Therefore this system applies to a whole variety 
of possible situations. 

We now turn our attention to the energy balance, that is we ask whether the above current- 
coupling system possesses a Hamiltonian structure. Remarkably, a positive answer is provided 
by the following Poisson bracket 

rr.^^ f \^F' 6G] . f ( 5F ^5G 5G 5F\ . 
{F,G}= m- -^,-^ d^x- /p — -V— - — -V— d^x 
J [dm omj J \oin op dm op J 

r f5F_ 5G__5F_ d_6G_ 6G_ _d_SF\ . 3 

J \ (5m 5m 6m dp 6f 6m dp 6f J 

f J{6F 5G\ „ d 6F d 5G\ ^. ,, 

B.(f xVx§-|^xVx§)d»x, (56) 

\om. 6d dm 6d J 

together with the Hamiltonian 

i/(m,p,/,B) = i [\j^d?^+±- //|p|2d3xd3p+ fpU{p)d^^ + ^ /iBl^d^x, 
2 y p 2mft J J 2po J 

(57) 

where the fluid velocity is replaced by the fluid momentum m = pu. The nature of the 
Hamiltonian structure (56) will be discussed later in relation to the MHD bracket (see theorem 
2). 

3.3 Pressure-coupling hybrid MHD schemes 

In this section we show how the Hamiltonian structure of the previous current-coupling scheme 
represents a helpful tool for the formulation of a pressure-coupling scheme. This scheme 
establishes an equation for the total momentum 

M = m -F / p / d^p , 
13 



under the assumption that the averaged kinetic momentum K = jp f d^p does not contribute 
to the total energy of the system. The definition of the momentum coordinate M yields the 
simple functional relations 

6F _ 6F 6F _ dF 5F _ 6F 

'6^~'5M' ^~^' Jb'Jb' 

which can be easily substituted in the bracket (56) to produce the new Hamiltonian structure 



6F _ 6F 6F 

If ~ 7/ ^^' 5M 



{F,G} 



M 



6F_ 5G 



d^x 



/ 6F 6G 



5G ^6F\ o 

— - • V^ d^x 
5M 6pJ 



+ 



+ 



f 
f 
B 



SF SG] ^ d 6F d 6G\ . . 

SF 6G\ (SG SF\\. ,3 



/ 6F ^ SG SG _ SF\ 
X V X X V X — 

V^M s^ sm s^j 



d^x 



(58) 



Let us now express the Hamiltonian (57) in terms of the total momentum M: we obtain 



2 



IM-KI 



■d-^x + 



1 



P 2m/,, 

Then, the Poisson bracket (58) yields the velocity equation 



/|p|2d3xd3p+ I pU{p)d^x+— I iBpd^x. (59) 

2/io 



— -+ w-V u = --Vp 
at p 



1 /■ o 1 

div / pp / dp H curl B x B , 

m/iP J Pop 



while the other equations remain identical to those in (53)-(55). Before proceeding further, we 
remark that neglecting all terms involving the averaged kinetic momentum Jp / d'^p (so that 
u\/u ~ U-'VU) produces the hybrid MHD model in [7] (although the general Vlasov equation 
is adopted here, rather than a gyrokinetic equation). However, this crucial step breaks the 
energy-conserving nature of the system and, when an energy balance equation is required, the 
model needs substantial modifications. 

3.3.1 First pressure-coupling MHD scheme 

A structure-preserving approximation can be easily obtained by neglecting the averaged kinetic 
momentum directly in the expression of the hybrid Hamiltonian (59) (so that M ~ m), which 
then becomes 



H{M,pJ,B) 



1 



- — ^d'^x + 



2m; 



|2j3 



/|prd^xd^p+ pU{p)d-^x + 



1 

2po 



iBpd^x 



Then, the Poisson bracket (58) yields the final equations 



f.,^^v,^ 



11/" 1 

--Vp div / pp/d^p H curlB x B 

P rr\hp J Pop 



-^+[/ + -^--^--^-V[/-p + a,,pxB--^ = 
at \ mhj ox ap op 

-J^ + div{pU) = , — = curl([/xB), 



(60) 

(61) 
(62) 
(63) 



where the term {/, p-C/} again appears to balance the pressure term in the energy conservation. 
Except for the Vlasov equation (62), the above system is identical to the scheme presented in 
[11, 36]. 
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3.3.2 Second pressure-coupling MHD scheme 

Notice that the following expression of the Haniiltonian is also possible 



H{M,p,f,B) 



|M| 



-d^x + 



1 



2m, 



/ 



K 

P 

n 



d^xd^p+ fpU{p)d^x+^ /jBj^d^x, 
J 2/io J 



(64) 



which in turn yields the equations 



dt 


- {u ■ v;c 




T'- 


u 

•^hp 


^V 


df 

dt ^ 


(P-V 
V ruh 


^ 
+ u 


5/ 

5x 


df / 
dp \ 


'du 

,9x 


dt ^ 


dw{pU) 


= 0, 


aB 

dt 


= cur 


(f/ 



dv p-v\ , ^,, „ a/ 

P + 7^ • ^^ + a/, p - y X B • -^ 

ax mh ) dp 



(65) 

(66) 

(67) 



where V = K/n. Notice that although complicated, the above kinetic equation yields the 
simple advection equation dfU + div(nt/) = for the zero-th moment n = Jfd^p. Then, 
the total density D = p + ruh n also satisfies dtD + div(D U) = 0, which is the total density 
equation appearing Cheng's hybrid model [7]. This transport equation for D and the velocity 
equation above both appear exactly the same in the hybrid scheme presented in [21, 42]. 
As a conclusion to this section, we remark that in the present Hamiltonian setting a purely 
advection equation for D is necessarily accompanied by the presence of the relative pressure 
tensor j {p — V) f d'^p in the velocity equation, instead of the absolute pressure / p® ^/ d^p. 
This fact is suggestive that the scheme in [21, 42] is preferable to the one presented in [7] 

3.4 The Hamiltonian structure and its consequences 

In the previous sections, the current-coupling hybrid kinetic scheme (52)-(55) was derived by 
inserting neutrality and Ohm's ideal law in the general Hamilton's equations arising from the 
Poisson bracket (45). Then, upon assuming negligible density and momentum for the energetic 
component, a simple momentum shift produced the pressure-coupling hybrid schemes. This 
section explains the origin of the hybrid Hamiltonian structures, starting from a direct sum 
Poisson bracket. In particular, the reason why the hybrid Poisson brackets (56) and (58) both 
satisfy the Jacobi identity is not at all evident and it needs particular care. 

We shall proceed by first deriving the bracket (56) from the direct sum bracket 



{F,G} 



N 



6F_ 5G_ 
5F 



d'^x 



P 
5G 



5F_ 6G 

In' J^ 



6F 



^ div -— - -I- -— - div ■ ^ 
(5m 5A (5N 6A 



Vx 



5G_ 5F_ 
6F_ 5G_ 

W^ Ja 



d'^x 



5G 5F 



d^x 



+ 



^(6F 5G\ n . 
f < -^,^ > d'^xd'^p. 
[6/ 6f 



(68) 



where the first two lines coincide with the bracket of ideal compressible MHD [17, 34, 31], 
while the last terin determines the Poisson bracket of the Vlasov equation. The momentum 
variable N is a new variable that will be given physical meaning in theorem 2. Notice that in 
the above bracket, the Vlasov distribution /(x, p) is expressed as a function of the mechanical 
momentum variable p = ruhV + QaA and not as a function of the variable p = ruhv. Thus one 
has 

/(x, p) = /(x, p + qhA) = /(x, p) 
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(Notice the different notation: tlie nieclianical momentum is denoted by the italic bold font p, 
while the roman bold notation p corresponds to nifiV.) 

Remark 6 (Lie algebra of the current-coupling scheme) The Poisson bracket (68) is a 
Lie-Poisson bracket on the direct-sum Lie algebra 

C°°{M.^) © (X{W^) (D (C°°(m3) 02(m3)) 

where C°°(M^) is the algebra of Hamiltonian functions (dual to Vlasov distributions), while the 
variables (N, p, A) belong to the dual space of the semidirect-product Lie algebra 

x{w^) (s)(c°°(]r3) © o2(]r3)) . 

Here O^(M^) is the space of differential two- forms containing the electric induction D • dS 
(where dS = h dS denotes the product of an infinitesimal surface element dS with its normal 
unit vector h). 

The Poisson bracket (56) is equivalent to (68) in the following sense: 

Theorem 2 The Poisson bracket (68) transforms to the bracket (56) under the momentum 
shift 

m = N + g;,A ^/(x, p) d^p = N + q^A J /(x, p) d^p . 

This result is proved in the appendix A.l. Perhaps the most direct physical consequence resides 
in the following circulation theorem 

Corollary 3 (Kelvin's theorem for the current-coupling scheme) For every loop 'yt{u) 
moving with velocity u = m/p, the current- coupling scheme (52) -(55) possesses the following 
circulation theorem 



d 



f (it-g/i — A)-dx=-® -Bx(/ig^VxB-a/i(K — mhuu)) ■ dx 

Jytiu) \ P J J'ytiu) P ^ ^ 



-\- Oh f - (V • (K — m/j nu) ) A • dx , 

J-rtiu) P 



where K = jp f d^p. 

Proof. From theorem 2, one concludes that the current-coupling scheme (52)-(55) possesses 
the Poisson bracket (68) and the Hamiltonian (57) in the form 

N + qhAJ fd^p 



27 p 



d^x + -^ I f \p- Qh A|2 d^xd^p 



2mh 



I p U{p) d^x + -^ / |V X Ap d^x . (69) 
J 2/xo J 



Then, replacing the Hamiltonian (69) in the bracket (68) yields the momentum equation 

- + /.JN=-pV- + -x(VxA)-AV.- (70) 

where £u denotes the Lie derivative along the velocity vector field u. After dividing by p and 
calculating 

— =qh -(N + qhnA)-^ p^'^V x V x A - a/j f{p-qh.A) d^p 

=Pq V x B — Oft (K — nih nu) , 

the proof follows easily by taking the circulation integral on both sides of (70). ■ 
Another relevant consequence of theorem 2 is the existence of Casimir functionals 
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Corollary 4 (Casimir functionals) The Poisson bracket (56) possesses the following Casimir 
functional 

c(/)=y<&(/)d3xd3p. 

Proof. Although it can be proven by direct verification, this is a direct consequence of the 
Casimir functional C{f) = f A{f) d^xd^p for the Poisson bracket (68), which is inherited from 
the general properties of the Hamiltonian structure of the Vlasov equation. ■ 

Notice that the magnetic helicity Ti = J A-V x A d^x is also a Casimir for the bracket (68), 
although it is not strictly a Casimir for the bracket (56), which is rather expressed in terms of 
the magnetic induction B = V x A. 

Analogous properties as those above also hold for the pressure coupling scheme. This can 
be seen by the following relation between the Poisson brackets (58) and (68): 

Theorem 5 The Poisson bracket (68) transforms to the bracket (58) under the change of 
variable 

M=N + Jp /(x, p) d^p = N + y p /(x, p) d^p + QhA I /(x, p) d^p . 

The proof of the above theorem is presented in appendix A. 2. Again, the above statement 
yields the 

Corollary 6 (Kelvin's theorems for pressure-coupling schemes) The Kelvin circulation 
laws for the pressure coupling schemes (61)-(63) and (65)-(67) read respectively as 

A / (u-- [pf{x,p) d^p") . dx = - / -Bx (^-ly X B - ahK) ■ dx 

+ ah(p - (V • K) A • dx , 

J-ytiU) P 



i f [/-- /p/(x,p)d^pj •dx = -/Zq^M p~^BxVxB-dx, 



di J^,t{u) \ P J J J-,t{U) 

where K = /p/d^p and the loop 'yt{U) moves with the velocity U = M/p 



Proof. Because of theorem 5, the pressure coupling scheme (61)-(63) can be written in terms 
of the momentum variable N = M — J'p/(x,p) d^p. Indeed, the equation for N is given 
immediately by the Poisson bracket (68) together with the Hamiltonian (60) in the form 



.-ifl 



2 



f-N+fp /(x, p) d^p d^x +-}- ff\p-q,^Af d^x d^p 
J P J 2mft J 

+ Ip U{p) d^x + -A / [V X Ap d=^x . (71) 
J ^Po J 

At this point, the Poisson bracket (68) produces the momentum equation (70), while the 
Hamiltonian (71) yields 

5H , 

—- = fio WxB-a/jK. 
OA 

Then, after replacing the above variational derivative in equation (70), the remaining steps 
coincide with those already followed in the proof of corollary 3. 
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The same arguments also hold for the second pressm'e-coupling scheme (65)-(67). In this case, 
the Hamiltonian (64) is written explicitly in terms of N as 



H 



1 r 1 



N + yp/(x,p)d^p 
1 r 1 



d''x + 



1 



2mA 



2j3, 



f\p-qhA\'d'^d'p 



2mh 



hJ Jf{-i^,p)^^p 



P /(X, p) d^p - QhA / /(x, p) d^p 

1 



+ pU{p)d-'x + 



2fio 



VxApd^x, (72) 



which in turn yields the following variational derivative: 



5H_ 

Ia 



-r^ =^0"! V X B. 



Then, inserting the following relation in the equation (70) arising from the Poisson bracket 
(68) and following the same steps as in the previous cases completes the derivation of the 
circulation law for the second pressure-coupling scheme (65)- (67). ■ 

Notice that if the hot particle density n is negligible at all times, then the above circulation law 
for the first pressure-coupling scheme (61)-(63) reduces to the circulation law given in corollary 
3 for the current coupling scheme. However, this particular case presents some obstacles that 
will be discussed in the next section. 

It is easy to verify that functionals of the type J'$(/)d^xd^p are also Casimirs of the 
Poisson bracket (58). The identification of Casimir functionals is a fundamental feature of the 
Hamiltonian approach, which allows for a nonlinear stability analysis as extensively described 
in [18]. The study of the stability properties of the hybrid MHD schemes is a considerable part 
of future work plans. 

Thus, we have the following consequence 

Corollary 7 The brackets (56) and (58) are Poisson brackets, whose Jacobi identity is inher- 
ited from the direct sum Poisson bracket 



Remark 7 (Lie algebra of pressure-coupling schemes) From an algebraic point of view, 
the Poisson bracket (58) is a Lie Poisson bracket [16] on the semidirect-product Lie algebra 

X(m3)(s)(C°°(]R3) X 01(]R3) X C7~(m6)) 



where we have denoted by X(M?) the space of velocity vector fields, while r2^(M'^) denotes the 
space of differential one-forms (i.e. the space of magnetic vector potentials). Notice that 
the space C°°(M^) is a Poisson algebra, which is dual to the Poisson manifold of Vlasov dis- 
tributions. This construction containing the semidirect-product Lie algebra X(M^)(DC°°(M^) 
appears in a continuum model for the first time. 

3.5 Arguments in favor of the current-coupling scheme 

The current coupling scheme for hybrid MHD is often considered to be equivalent to ordinary 
pressure coupling schemes [36, 35, 11]. However, the assumption of a rarefied energetic ion 
plasma, poses nontrivial problems in the mathematical formulation. Indeed, we recognize that 
the assumption of negligible particle density n may lead to some ambiguities. For example, we 
saw that in the second pressure-coupling scheme (65)-(67), the condition n = is preserved 
by the dynamics, although this reflects in a possible divergence of the Hamiltonian, which 
indeed contains the term n~ |K| . Then, the only possibility to avoid this divergence is to set 
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K = 0. For example, the whole velocity K/n is sometimes neglected completely [42]. However, 
the well known pressure term divj" (p — V) /d^p appearing in the equation for K forbids 
this possibility unless one accepts that this pressure term is also negligible, which contradicts 
the starting hypothesis of non-negligible pressure effects in the total momentum dynamics. 
This non-negligible pressure term is a continuous source of hot particle averaged momentum, 
whose initial low levels are not maintained by the dynamics. Notice that this phenomenon is 
not peculiar of the pressure-coupling schemes presented in this paper; rather, this problem is 
intrinsically unavoidable and it is present in all pressure-coupling schemes. 

The current-coupling scheme provides an exceptional way of avoiding this sort of problems, 
since it makes no assumption on the moments of the Vlasov distribution. Some advantages 
of the current-coupling scheme were also emphasized in [1]. Moreover, the existence of the 
newly discovered Hamiltonian structure for this model, together with its consequent circulation 
theorem and Casimir functions, enriches this model significantly. In conclusion, the pressure 
coupling scheme is very different from the current coupling scheme and the latter solves all the 
significant problems emerging from pressure coupling schemes. 

4 Hybrid kinetic-fluid formulation of Hall MHD 

In the preceding sections we have investigated the hybrid kinetic-fluid formulation of MHD 
models taking into account the dynamics of energetic ions. On the other hand, the existence of 
energetic ions is also a peculiar feature of Hall effects in magnetized plasmas. Indeed, one of the 
main consequences of Hall effects is the progressive acceleration of ions that become particularly 
fast. In Hall MHD models, this scenario leads to the simple idea of treating the whole ion 
dynamics by a kinetic approach, while the electrons are still considered as a neutralizing 
background. This model is usually referred to as "kinetic-particle-ion/fluid electron" hybrid 
scheme [35, 45, 4, 47]. Its equations are easily derived from the kinetic- multifluid system of 
section 3.1, upon letting the only label s = 2 identify the electron species. Then, the plasma 
equations of motion reduce to the Vlasov kinetic equation (41) and the electron fluid equation 
(39) (with s = 2). Neglecting electron inertia (that is letting me -^ 0) and assuming neutrality 
(that is letting eo — > 0) yields the well known equations [35, 45, 4] 



df V df ^ -D ^/ , 9i 



dt rrii 5x 9p aePe 



fleVpe - iaipf d^p - /Xq ^ V X B j X B 



-^ B X U, /p / d^p - ^-1 V X B j 



a/ _ 

dp 

(73) 

(74) 



In the search for a Hamiltonian structure of the above equations, it is helpful to recall the 
Hamiltonian structure of Hall MHD, which was presented in [IG]. Proceeding by analogy to 
the Hall-MHD bracket in [IG], one introduces the Poisson bracket 

where Qe = UePe is the electron charge density and the kinetic Vlasov term in the last line 
replaces the fluid bracket terms appearing in the Hamiltonian structure of Hall MHD [16]. 
Notice that the first line of (75) is equivalent to a fluid bracket in the density-momentum 
variables (Qe,QeA) [IG]. At this point, the Hamiltonian structure of the equations (73)-(74) 
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arises by simply expressing the above bracket in terms of the Vlasov distribution /(x, p) = 
/(x,p — QiA) = f{x,p). Then, the final result is 

,^^, /"/^ 1^ ,^ SF 6G . [(5F ^5G 5G ^5F\n 

{F,G}=y(Q-VxA).-x-d3x-y(^-.V— --.V— jd3x 

[J{5F5G\ ^ , d 5F d 6G\ . . 

which is indeed the Poisson bracket of equations (73)-(74) and it is accompanied by the Hamil- 
tonian 

if(/,A,Qe) = -^ [f\p\'d'Kd'p+ L(Q,)d3x+-^ /iVxApd^x, (77) 

^rrij J J ZfiQ J 

where (piQe) appropriately denotes the internal energy expressed in terms of the charge Qe- 
Notice that the charge neutrality Qe + qiU = is preserved by the dynamics, as it is shown by 
the equation 

dQe /" s 9 f r, 

-^ = aidiv p f d p = -Qi — p f d p , 



which arises from the Hamiltonian structure (76). 

It is also interesting to notice that the class of Casimir functionals C(/) = /A(/) d^xd^p 
for the Vlasov equation again produces the Casimir functionals C(/) = f <!*(/) d'^xd'^p for the 
Hamiltonian structure (75) under consideration. Here, the magnetic helicity 

■H= / A- V X Ad^x 



is a Casimir of both Poisson brackets (75) and (76). These Casimir functionals provide the 
opportunity for a nonlinear stability analysis, which will the the subject of future work. 

Remark 8 (Lie algebraic structure) The Hamiltonian structure (76) of the hybrid Hall- 
MHD model is inherited from the Lie-Poisson bracket (75) on the direct-sum Lie algebra 

'^^®G°°{R^)) eC7~(M^), 



where C°°{M.) is the Poisson algebra of Hamiltonian functions, while the variables (QgA, Qe) 
belong to the dual space of the semidirect-product Lie algebra X(M^) (S)C°°(M^) [16]. The 
Hamiltonian corresponding to the Poisson bracket (75) is given by 

H{f,A,Qe) = ^ [f\p-qiAfd^xd^p+ /'^(Q,) d^x + -^ /iVxApd^x. (78) 

The absence of a fluid velocity equation in the hybrid formulation of Hall MHD makes this 
model particularly difficult to implement numerically. However, one can introduce a velocity 
equation by simply coupling the ion kinetic equation to the equation for its first two moments 
{piVi, Pi) = (Jp / d^p, rrij // d^p). Evidently, this operation does not produce any changes to 
the physics, as long as the form of the Hamiltonian is not modified. On the other hand, the 
ion momentum equation might be useful to isolate kinetic pressure terms explicitly [35]. When 
these kinetic pressure terms are omitted, the ion momentum equation consistently returns the 
ordinary equations of Hall MHD. 
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5 Conclusions and future work plans 

This paper has presented the Hamiltonian structure of several hybrid kinetic-fluid models. 
Upon making assumptions on the form of the Hamiltonian, new pressure-coupling schemes 
were derived for either ordinary fluids and MHD, whereas the existing current-coupling MHD 
scheme [36, 35] was shown to possess an intrinsic Hamiltonian structure. All Poisson brackets 
were derived from first principles, so that the Jacobi identity is always guaranteed. In turn, the 
Hamiltonian structure produces a Casimir functional and a circulation law for each of these 
models. The same approach was used to present the Hamiltonian structure of the existing 
hybrid Hall-MHD scheme [45, 47]. 

The next steps will focus on the nonlinear stability analysis, which is now made possible 
by the newly discovered Hamiltonian structures. Another relevant open question concerns 
the extension of these Hamiltonian structures to the case of gyrokinetic equations, which are 
frequently used in simulations. The Hamiltonian structure of gyrokinetic equations is well 
known [2, 23, 13] and it is possible that the same methods used in this paper can be fruitful 
also in that context. 

While the Hamiltonian approach was used in this paper, its Lagrangian counterpart may 
also deserve further investigation. This direction can be easily approached by the methods 
developed in [5] and is part of current work [20] . 

As a conclusive remark, it is worth mentioning that while pressure coupling schemes were 
shown to require transport-like terms in the accompanying kinetic equation, one cannot exclude 
a priori that removing these terms still allows for other exotic Poisson brackets. Indeed, 
although this is unlikely, it cannot be excluded. However, even if this is the case, the alternative 
Poisson brackets would not be of Lie-Poisson form and thus they would not possess well defined 
symmetry properties, which are the main focus of this paper. 
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A Appendix 

A.l Proof of theorem 2 



First, one computes the variational derivatives 
6F _ 6F 5F 

Then, the Poisson bracket (68) becomes 
' 5F_ 6G ' 

5F 



^ + ^'^fciy^("'^)^^' Jf=Jf 



6F ^ 6F 

om 



{F,G} 



N- 
A 



5G 



6F 8G 
om op 



5m 5p J 



'''\ s. 



dm oA 



5G , 5F 
om oA 



V 



5F_ 5G__5G_ 5F_ 
5vii 5A. 5vii 5A. 



d^x 



Qh 



5F ^ f 5G\ 
-— div n-— 
om \ Omy 



5G , / 5F\ 
-—d\Y[n—- 
Om \ omy 



2V X n 



5F_ 5G_ 
(5m (5m 



d'^x 



+ 



4f-f}-^'--'/^({f-a-{f-a)-^^'- 



where we have dropped the bar notation for convenience and we have denoted n = J f d^p. 
Upon expanding the third line above, one has 
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J div ■ 
5m oA 
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5m oA 



\ (3m 
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5F 
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Vn 
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d^xd^p 
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d^xd^p 



Then, the standard vector identities 

a X b X c = (a • c) b — (a • b) c 

V x (a X b) = a (V • b) - b (V • a) + [a, b] 



yield 
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Upon collecting all vector field commutators and simplifying according to the general formula 

Vx{ka) = kVxa + Vk x a , 
we obtain 
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m 
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6F 5G' 
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At this point, one expresses the above Poisson bracket in terms of the distribution /(x, p) 
/(x,p — QhA.). Then, the relations [25, 27] 



r f ," 1 r , , 1 ^ ^ dh dk 
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transform the bracket to 
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Introducing the magnetic induction B = V x A by the formula 



51^ 

Ja 



Vx 



5£ 



completes the proof, thereby yielding the bracket (56). 



A. 2 Proof of theorem 5 

The momentum shift M = N + Jp f d'^p takes the Poisson bracket (68) to the form 
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where £u denotes Lie derivative with respect to the velocity u. Then, substitution of the 
formulas 



{/i,fc} = {/i,A;} + <?^VxA-(|^x|^^ 



5F 5F (' d 6F . 



and some integration by parts yield the bracket 
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Then, using the formula 

yields the Poisson structure (58). 



6F _ 6F 



(83) 

(84) 
(85) 

(86) 
(87) 



(88) 



24 



References 

[1] Belova E.V.; Denton R.E.; Chan A. A. Hybrid simulations of the effects of energetic par- 
ticles on low-frequency MHD waves. J. Comput. Phys. 136 (1997), no. 2, 324-336 

[2] Brizard, A.J.; Hahm, T.S. Foundations of nonlinear gyrokinetic theory. Rev. Mod. Phys. 
79 (2007), no. 2, 421-468 

[3] Brizard, A.J.; Tracy, E.R. Mini- conference on Hamiltonian and Lagrangian methods in 
fluid and plasma physics. Phys. Plasmas 10 (2003), no. 5, 2163-2168 

[4] Byers, J. A.; Cohen, B.I.; Condit, W.C.; Hanson, J.D. Hybrid simulations of quasineutral 
phenomena in magnetized plasma. J. Comp. Phys. 27 (1978), no. 3, 363-396 

[5] Cendra, H.; Holm, D.D.; Hoyle, M.J.W.; Marsden, J.E. The Maxwell-Vlasov equations in 
Euler-Poincar form. J. Math. Phys. 39 (1998), no. 6, 3138-3157 

[6] Chandy, A.J.; Frankel, S.H. Regularization-based sub-grid scale (SGS) models for large 
eddy simulations (LES) of high-Re decaying isotropic turbulence. J. Turbul. 10 (2009), no. 
25, 1-22 

[7] Cheng, C.Z. A kinetic-magnetohydrodynamic model for low-frequency phenomena. J. Geo- 
phys. Res. 96 (1991), no. A12, 21,159-21,171 

[8] Cheng, C.Z.; Johnson, J.R. A kinetic-fluid model. J. Geophys. Res. 104 (1999), no. Al, 
413-427. 

[9] Fla, T. A hybrid fluid-kinetic theory for plasma physics. arXiv:math-ph/0110036vl 

[10] Foias, C; Holm, D.D.; Titi, E.S. The Navier- Stokes- alpha model of fluid turbulence. Phys. 
D 152/3 (2001) 505519 

[11] Fu, G.Y.; Park, W. Nonlinear hybrid simulation of the toroidicity-induced Alfven Eigen- 
mode. Phys. Rev. Lett. 74 (1995), no. 9, 1594-1596 

[12] Gibbons, J.; Holm, D.D.; Tronci, C. Geometry of Vlasov kinetic moments: a bosonic Fock 
space for the symmetric Schouten bracket. Phys. Lett. A 372 (2008), 4184-4196 

[13] Grebogi, C; Kaufman, A.N.; Littlejohn, R.G. Hamiltonian theory of ponderomotive effects 
of an electromagnetic wave in a nonuniform magnetic field. Phys. Rev. Lett. 43 (1979), 
no. 22, 1668-1671. 

[14] Holm, D.D. Hamiltonian dynamics of a charged fluid, including electro- and magnetohy- 
drodynamics. Phys. Lett. A 114 (1986), no. 3, 137-141. 

[15] Holm, D.D. Hamiltonian dynamics and stability analysis of neutral electromagnetic fluids 
with induction. Phys. D 25 (1987), no. 1-3, 261-287. 

[16] Holm, D.D. Hall magnetohydrodynamics, conservation laws and Lyapunov stability. Phys 
Fluids 30 (1987), no. 5, 1310-1322. 

[17] Holm, D.D.; Kupershmidt, B.A. Noncanonical Hamiltonian formulation of ideal magne- 
tohydrodynamics. Phys. D 7 (1983), no. 1-3, 330-333. 

[18] Holm, D.D.; Marsden, J.E.; Ratiu, T.S.; Weinstein, A. Nonlinear stability of fluid and 
plasma equilibria. Phys. Rep. 123 (1985), no. 1-2, 1-116 



25 



[19] Holm, D.D.; Tronci, C. Geodesic Vlasov equations and their integrable moment closures. 
J. Geom. Mech. 1 (2009), no. 2, 181- 

[2o; 

[21 



[22; 

[23 
[24 
[25; 
[26 
[27; 

[28 

[29 

[30 

[31 

[32; 
[33 

[34; 

[35; 

[36 

[37; 



Holm, D.D.; Tronci, C. Euler-Poincare formulation of hybrid plasma models. Preprint 

Kim, C.C.; Sovinec, C.R.; Parker, S.E.; the NIMROD team Hybrid kinetic-MHD simula- 
tions in general geometry. Comp. Phys. Comm. 164 (2004), 448-455 



Liboff, R.L.; Kinetic Theory, Springer-Verlag 1998. 

Littlejohn, R.G. Variational principles of guiding centre motion. J. Plasma Phys. 29 
(1983), no. 1, 111-125 

Low, F.E. A Lagrangian formulation of the Boltzmann- Vlasov equation for plasmas. Proc. 
R. Soc. London, Ser. A 248 (1958), 282-287 

Marsden, J.E.; Weinstein, A. The Hamiltonian structure of the Maxwell- Vlasov equations. 
Phys. D 4 (1981/82), no. 3, 394-406 

Marsden, J.E.; Weinstein, A.; Ratiu, T.S. Semidirect products and reduction in mechanics. 
Trans. Amer. Math. Soc. 281 (1984), no. 1, 147-177. 

Marsden, J.E.; Weinstein, A.; Ratiu, T.; Schimd, R.; Spencer, R.G. Hamiltonian systems 
with symmetry, coadjoint orbits and plasma physics, Atti Accad. Sci. Torino CI. Sci. Fis. 
Mat. Natur. 117 (1983), no. 1, 289-340 

Morrison, P.J. Hamiltonian and action principle formulations of plasma physics. Phys. 
Plasmas 12 (2005), 058102 

Morrison, P.J. Hamiltonian field description of the one- dimensional Poisson-Vlasov sys- 
tem. Princeton Plasma Physics Laboratory Report, PPPL-1788 (1981). 

Morrison, P.J. Hamiltonian field description of two-dimensional vortex fluids and guiding 
center plasmas. Princeton Plasma Physics Laboratory Report, PPPL-1783 (1981). 

Morrison, P.J. On Hamiltonian and action principle formulations of plasma dynamics. 
AIP Conf. Proc, 1188 (2009), 329-344. 

Morrison, P.J. Poisson brackets for fluids and plasmas. AIP Conf. Proc. 88 (1982), 13-46 

Morrison, P.J. The Maxwell- Vlasov equations as a continuous Hamiltonian system. Phys. 
Lett. A 80 (1986), no. 5-6, 383-386 

Morrison, P.J.; Greene, J.M. Noncanonical Hamiltonian density formulation of hydrody- 
namics and ideal magnetohydrodynamics. Phys. Rev. Lett. 45 (1980), 790-794. 

Park, W.; Belova, E.V.; Fu, G.Y.; Tang, X.Z.; Strauss, H.R.; Sugiyama, L.E. Plasma 
simulation studies using multilevel physics models. Phys. Plasmas 6 (1999), no. 6, 1796- 
1803. 

Park, W.; Parker, S.; Biglari, H.; Chance, M.; Chen, L.; Cheng, C.Z.; Hahm, T.S.; Lee, 
W.W.; Kulsrud, R.; Monticello, D.; Sugiyama, L.; White, R. Three-dimensional hybrid 
gyrokinetic-magnetohydrodynamics simulation. Phys. Fluids B 4 (1992), no. 7, 2033-2037 

Pfirsch, D. New variational formulation of Maxwell- Vlasov and guiding center theories 
local charge and energy conservation laws. Z. Naturforsch. A 39a (1984), 1-8 



26 



[38] Pfirsch, D.; Morrison, P.J. Local conservation laws for the Maxwell-Vlasov and collisionless 
kinetic guiding- center theories. Phys. Rev. A 3 (1985), no. 2, 1714-1721 

[39] Pfirsch, D.; Morrison, P.J. The energy-momentum tensor for the linearized MaxwellVlasov 
and kinetic guiding center theories. Phys. Fluids B 3 (1991), no. 2, 271-283 

[40] Spencer, R.G. The Hamiltonian structure of multi- species fluid electrodynamics. AIP Conf. 
Proc, 88 (1982), 121-126. 

[41] Spencer, R.G.; Kaufman, A.N. Hamiltonian structure of two-fluid plasma dynamics. Phys. 
Rev. A (3) 25 (1982), no. 4, 2437-2439 

[42] Takahashi, R.; Brennan, D. P.; Kim, C.C. Kinetic effects of energetic particles on resistive 
MHD stability. Phys. Rev. Lett. 102 (2009), 135001 

[43] Todo, Y.; Sato, T.; Watanabe, K.; Watanabe, T. H.; Horiuchi, R. Magnetohydrodynamic 
Vlasov simulation of the toroidal Alfven eigenmode. Phys. Plasmas 2 (1995), no. 7, 2711- 
2716. 

[44] Weinstein, A. Gauge groups and Poisson brackets for interacting particles and fields. AIP 
Conf. Proc, 88 (1982), 1-11. 

[45] Winske, D.; Yin, L.; Omidi, N.; Karimbadi, H.; Quest, K. Hybrid simulation codes: past, 
present and future - a tutorial. Lect. Notes Phys. 615 (2003), 136-165 

[46] Ye, H.; Morrison, P.J. Action principles for the Vlasov equation. Phys. Fluids B 4 (1992), 
no. 4, 771-777 

[47] Yin, L.; Winske, D.; Gary, S.P.; Birn, J. Hybrid and Hall-magnetohydrodynamics simu- 
lations of collisionless reconnection: effect of the ion pressure tensor and particle Hall- 
magnetohydrodynamics. Phys. Plasmas 9 (2002), no. 6, 2575-2584 



27 



